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Models of dark energy in which neutrinos interact with the scalar field supposed to be responsible 
for the acceleration of the universe usually imply a variation of the neutrino masses on cosmological 
time scales. In this work we propose a parameterization for the neutrino mass variation that captures 
the essentials of those scenarios and allows to constrain them in a model independent way, that is, 
without resorting to any particular scalar field model. Using WMAP 5yr data combined with the 
matter power spectrum of SDSS and 2dFGRS, the limit on the present value of the neutrino mass 
is mo = m v (z = 0) < 0.43 (0.28) eV at 95% C.L. for the case in which the neutrino mass was lighter 
(heavier) in the past, a result competitive with the ones imposed for standard (i.e., constant mass) 
neutrinos. Moreover, for the ratio of the mass variation of the neutrino mass Am„ over the current 
mass mo we found that log[| Am„ |/mo] < —1.3 (—2.7) at 95% C.L. for Am„ < (Am„ > 0), totally 
consistent with no mass variation. These stringent bounds on the mass variation are not related to 
the neutrino free-streaming history which may affect the matter power spectrum on small scales. 
On the contrary, they are imposed by the fact that any significant transfer of energy between the 
neutrino and dark energy components would lead to an instability contradicting CMB and large 
scale structure data on the largest observable scales. 

PACS numbers: 14.60.St, 98.80.-k, 98.80.Cq, 98.80.Es 



I. INTRODUCTION 

Since the accelerated expansion of the universe was 
first observed with Type la supernovae (SN) [l|, Q , the 
case for a cosmological constant-like fluid that dominates 
the energy density of the universe has become stronger 
and is well established by now with the new pieces of 
data gathered @ . 

Several candidates for the accelerating component of 
the universe, generically dubbed dark energy (DE), have 
been proposed [!, 0, 0, @ , but understanding them the- 
oretically and observationally has proven to be challeng- 
ing. On the theoretical side, explaining the small value of 
the observed dark energy density component, ~ (10 -3 
cV) 4 , as well as the fact that both dark energy and mat- 
ter densities contribute significantly to the energy bud- 
get of the present universe requires in general a strong 
fine tuning on the overall scale of the dark energy mod- 
els. In the case in which the dark energy is assumed 
to be a scalar field <f> slowly rolling down its flat po- 
tential V((f>), the so-called quintessence models Q, the 
effective mass of the field has to be taken of the order 
m = IcfVi^/dcj) 2 ] 1 / 2 ~ 10" 33 cV for fields with vac- 
uum expectation values of the order of the Planck mass. 

On the observational side, choosing among the dark 
energy models is a complicated task Q. Most of them 
can mimic a cosmological constant at late times (that is, 
an equation of state = V<t>l P<t> = — 1) @, an( i au data 
until now are perfectly consistent with this limit. In this 
sense, looking for different imprints that could favor the 
existence of a particular model of dark energy is a path 



worth taking. 

Our goal in this paper consists in understanding 
whether the so-called Mass- Varying Neutrinos (MaVaNs) 
scenario [l(| [HI, El, G3 could be constrained not only 
via the dark energy effects, but also by indirect signs of 
the neutrino mass variation during cosmological evolu- 
tion, since neutrinos play a key role in several epochs 
[ill, [IH . An indication of the variation of the neutrino 
mass would certainly tend to favor this models (at least 
on a theoretical basis) with respect to most DE mod- 
els. One should keep in mind that MaVaNs scenarios 
can suffer from stability issues for the neutrino perturba- 
tions [ljj, although there is a wide class of models and 
couplings that avoid this problem [H, [H, US [HI, H3 ■ 

Similar analyses have been made in the past, but they 
have either assumed particular models for the interac- 
tion between the neutrinos and the DE field [H, HH HH , 
or chosen a parameterization that does not reflect the 
richness of the possible behavior of the neutrino mass 
variations [26]. 

In order to be able to deal with a large number of 
models, instead of focusing on a particular model for the 
coupling between the DE field and the neutrino sector, 
we choose to parameterize the neutrino mass variation 
to place general and robust constraints on the MaVaNs 
scenario. In this sense, our work complements previous 
analyses by assuming a realistic and generic parameteri- 
zation for the neutrino mass, designed in such a way to 
probe almost all the different regimes and models within 
the same framework. In particular, our parameterization 
allows for fast and slow mass transitions between two 
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values of the neutrino mass, and it takes into account 
that the neutrino mass variation should start when the 
coupled neutrinos change their behavior from relativis- 
tic to nonrelativistic species. We can mimic different 
neutrino-dark energy couplings and allow for almost any 
monotonic behavior in the neutrino mass, placing reliable 
constraints on this scenario in a model independent way. 

Our work is organized as follows: in Scction|TT]we give a 
brief review of the MaVaNs scenario and its main equa- 
tions. In Section IIIII we present our parameterization 
with the results for the background and the perturbation 
equations obtained within this context. The results of 
our comparison of the numerical results with the data 
and the discussion of its main implications are shown in 
section IIVI Finally in section [V] the main conclusions 
and possible future directions are discussed. 



II. MASS- VARYING NEUTRINOS 



Since the total energy momentum tensor is conserved, 
the dark energy fluid equation also presents an extra 
right-hand side term proportional to the neutrino energy 
momentum tensor trace, = (p v — 3p„), and can be 

written as 



P4, + 'iHp,), (1 + w$) = -a(4>)(j) {p u - 3jj„) 



(5) 



For a homogeneous and isotropic scalar field, the energy 
density and pressure are given by 



2a 2 



V{<j>) 



P* = ^2- V (<f>), (6) 



and both equations lead to the standard cosmological 
Klein-Gordon equation for an interacting scalar field, 
namely. 



= -a 2 a(0) {p v - 3p„) . (7) 



In what follows, we consider a homogeneous and 
isotropic universe with a Robertson- Walker flat metric, 
ds 2 = a 2 (dr 2 + dr 2 + r 2 dfi 2 ) , where r is the confor- 
mal time, that can be written in terms of the cosmic 
time t and scale factor a as dr = dt/a, in natural units 
(H = c = fcs = 1). In this case, the Friedmann equations 
read 



n 1 



n = - 



6m 2 



3m 2 



(P + 3p) ; 



(1) 

(2) 



where the dot denotes a derivative with respect to con- 
formal time, and the reduced Planck mass is m p = 
1/V8^G = 2.436 x 10 18 GeV. As usual, p and p cor- 
respond to the total energy density and pressure of the 
cosmic fluid, respectively. The neutrino mass in the mod- 
els wc arc interested in is a function of the scalar field (f> 
that plays the role of the dark energy and can be written 
as 



m„(0) =M„/(0) 



(3) 



where M v is a constant and different models are repre- 
sented by distinct /(</>). 

The fluid equation of the neutrino species can be di- 
rectly obtained from the Boltzmann equation for its dis- 
tribution function I24II. 



p v + 3Hp v (1 + w v ) = a{4>)(j) (p v - 3p u ) 



(4) 



where oc((j>) = dln[m v ((j))]/d(j) takes into account the vari- 
ation of the neutrino mass, and w x = p x j 1 p x is the equa- 
tion of state of the species x. For completeness and later 
use, we will define Q, x q = p x /p c o, the standard density 
parameter, where the current critical density is given by 
Pen = 3H^m 2 p = 8.099 h 2 x 1CT 11 eV 4 and H = 100 h 
km s _1 Mpc -1 is the Hubble constant. 



From the above equations one sees that, given a potential 
V((j)) for the scalar field and a field-dependent mass term 
m u {4>) for the neutrino mass, the coupled system given 
by equations (Q]), (g]), and ([7]), together with the fluid 
equations for the baryonic matter, cold dark matter and 
radiation (photons and other massless species) can be 
numerically solved [24|. Notice that a similar approach 
has been used for a possible variation of the dark matter 
mass 1271 and its possible interaction with the dark en- 
ergy [28l. [29|. with several interesting phenomenological 
ramification s [3Cl.[3ll . l32l . l33l . HH, [3f| . 

Following [3l|, [33| , equations ((U) and (O can be rewrit- 
ten in the standard form. 



Pv 



3U Pv 1 + < 



(off; 



37Y/90 1 + u> 



„(eff) 



= , 
= , 



(8) 



if one defines the effective equation of state of neutrinos 
and DE as 



■w 



(off) 



„(off) 



Pv_ _a_ 

Pv 

P<t> . 



{pu - 3p v ) 



3Hp v 
{p v - 3p v ) 



(9) 



3Hp<f> 



The effective equation of state can be understood in 
terms of the dilution of the energy density of the species. 
In the standard noncoupled case, the energy density of a 
fluid with a given constant equation of state w scales as 
p oc a~ 3 ( 1+ "'). However, in the case of interacting fluids, 
one should also take into account the energy transfer be- 
tween them, and the energy density in this case will be 
given by 



p(z) = po exp 



l + w (0ff) (2')Wi(l + 2') 



(10) 



where the index denotes the current value of a parame- 
ter, and the redshift z is defined by the expansion of the 
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scale factor, a = ao(l + z) 1 (in the rest of this work we 
will assume a = 1). For a constant effective equation of 

state one obtains the standard result, p oc a~ 3 ( 1+v)l \ 
as expected. 

Notice that this mismatch between the effective and 
standard DE equations of state could be responsible for 
the "phantom behavior" suggested by supernovae data 
when fitting it using a cosmological model with nonin- 
teracting components (3^|. This effect could be observ- 
able if dark energy was coupled to the dominant dark 
matter component. For the models discussed here, how- 
ever, it cannot be significant: the neutrino fraction today 
(U^o/n^o ~ 10~ 2 ) is too small to induce an "effective 
phantom-like" behavior. 

As we commented before, the analysis until now dealt 
mainly with particular models, that is, with particular 
functional forms of the dark energy potential V(<j>) and 
field dependence of the neutrino mass a(4>). A notice- 
able exception is the analysis of Ref. [2(|, in which the 
authors use a parameterization for the neutrino mass a 
Id Chevallier-Polarski-Linder (CPL) [1, [H, [H : m v (a) = 
m v o+m v \(l— a). However, although the CPL parameter- 
ization works well for the dark energy equation of state, 
it cannot reproduce the main features of the mass varia- 
tion in the case of variable mass particle models. In the 
case of the models discussed here, for instance, the mass 
variation is related to the relativistic/nonrelativistic na- 
ture of the coupled neutrino species. With a CPL mass 
parameterization, the transition from mi to too always 
takes place around z ~ 1, which is in fact only compat- 
ible with masses as small as 10 -3 eV. Hence, the CPL 
mass parameterization is not suited for a self-consistent 
exploration of all interesting possibilites. 

One of the goals in this paper is to propose and test 
a parameterization that allows for a realistic simulation 
of mass-varying scenarios in a model independent way, 
with the minimum possible number of parameters, as ex- 
plained in the next section. 



III. MODEL INDEPENDENT APPROACH 



A. Background equations 



momentum, but not on the mass [161 ]. In what follows 
we have neglected the small spectral distortions arising 
from non- instantaneous neutrino decoupling [38|]. Thus, 
the neutrino energy density and pressure are given by 



1 



Pv 



dq 



dnq 2 ef°(q) , 



(12) 



V» = h\ iXdnqifiqf-, (13) 
3a 4 J {2tt) ,s e 

where e 2 = q 2 + m 2 (a)a 2 (assuming that m u depends 
only on the scale factor). Taking the time-derivative of 
the energy density, one can then obtain the fluid equation 
for the neutrinos, 

p v + m (p v + p v ) = dln ™^ n {pv _ 3^) , (14) 

du 

where u = In a = — ln(l + z) is the number of e-folds 
counted back from today. Due to the conservation of 
the total energy momentum tensor, the dark energy fluid 
equation is then given by 

<ilnTO„(u) 

P0 + 3H/90 (1 + -u^) = — Tl{p u — 3p v ) . (15) 

We can write the effective equations of state, defined in 
eqs. ©, as 



w 



off 



p v c?lnTO„(u) / 1 
3 



Pv 



_ Pv_ 

du \ 3 p v 

d\a.m v (u) (\ p v 
ft(i> J du \3 p v 



(16) 



The above results only assume that the neutrino mass 
depends on the scale factor a, and up to this point, we 
have not chosen any particular parameterization. Con- 
cerning the particle physics models, it is important to 
notice that starting from a value of w^, and a function 
TOy(a) one could, at least in principle, reconstruct the 
scalar potential and the scalar interaction with neutrinos 
following an approach similar to the one in Ref. [3^ . 



As usual, the neutrino energy density and pressure are 
given in terms of the zero order Fermi-Dirac distribution 
function by 



f(l) 



<),■ 



(11) 



where q = ap denotes the modulus of the comoving mo- 
mentum qt = qm (S^mrij = 1), g v corresponds to the 
number of neutrino degrees of freedom, and T u q is the 
present neutrino background temperature. Notice that 
in the neutrino distribution function we have used the 
fact that the neutrinos decouple very early in the history 
of the universe while they are relativistic, and therefore 
their equilibrium distribution depends on the comoving 



B. Mass variation parameters 

Some of the main features of the MaVaNs scenario are: 
(i) that the dark energy field gets kicked and moves away 
from its minimum (if to^ > H ) or from its previous slow- 
rolling trajectory (if to^ < H) when the neutrinos be- 
come non-relativistic, very much like the case when it is 
coupled to the full matter content of the universe in the 
so-called chameleon scenarios [4(j ; and (ii) that as a con- 
sequence, the coupling with the scalar field generates a 
neutrino mass variation at that time. Any parameteriza- 
tion that intends to mimic scalar field models interacting 
with a mass- varying particle (neutrinos, in our case) for 
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the large redshift range to which the data is sensitive 
should at least take into account those characteristics. 
Moreover, the variation of the mass in most models (see 
j24| . for instance) can be well approximated by a tran- 
sition between two periods: an earlier one, in which the 
mass is given by mi , and the present epoch, in which the 
mass is given by mo (we will not consider here models in 
which the neutrino mass behavior is nonmonotonic) . The 
transition for this parameterization, as mentioned before, 
starts when neutrinos become nonrelativistic, which cor- 
responds approximately to 



1.40 



1 eV 



3 T, 



70 



mi 
1 cV 



2 x lrr 



mi 
FeV 



(17) 



where mi corresponds to the mass of the neutrino dur- 
ing the period in which it is a relativistic species. Before 
znr we can treat the neutrino mass as essentially con- 
stant, since the right-hand side (RHS) of the fluid equa- 
tion is negligible compared to the left-hand side (LHS), 
and therefore there is no observable signature of a possi- 
ble mass variation. 

When the neutrinos become nonrelativistic, the RHS 
of the DE and neutrino fluid equations becomes impor- 
tant, and the neutrino mass starts varying. In order to 
model this variation, we use two parameters, namely the 
current neutrino mass, mo, and A, a quantity related to 
the amount of time that it takes to complete the tran- 
sition from mi to mo. That behavior resembles very 
much the parameterization of the dark energy equation 
of state discussed in [4l[ , except for the fact that in our 
case the transition for the mass can be very slow, tak- 
ing several e-folds to complete, and must be triggered by 
the time of the nonrelativistic transition, given by equa- 
tion {T7J). Defining / = [1 + e"[" U+A)-u NR ]/A]-i and 
/* = [1 + e UNR ' A ] _1 we can use the form 



m + (mi - m ) x T(u, u NR , A) 



where 

r(u,ii NR , A) 



f_ 

u 



1 - 



1 



u N r/A 



1 + e -[« (l+A)-u NR ]/A 



(18) 



(19) 



Starting at unr = — ln(l + znr), the function 
r(it, unr, A) decreases from 1 to 0, with a velocity that 
depends on A. The top panel in Figure Q] gives the be- 
havior of eq. (fl~8|) with different parameters; the bottom 
panels shows that in this parametrization, the derivative 
of the mass with respect to c-fold number resembles a 
Gaussian function. The peak of the quantity dm/du oc- 
curs at the value u = wnr/(1 + A); hence, for A « 1. 
the mass variation takes place immediately after the non- 
relativistic transition (u ~ unr) and lasts a fraction of 
c-folds (roughly, 3A e-folds); for 1 < A < |wnr| the 
variation is smooth and centered on some intermediate 
redshift between znr and 0; while for A 3> |ujvi?|, the 




FIG. 1: (Color online) Neutrino mass behavior for the param- 
eterization given by equation (|18[) . Top panel: Neutrino mass 
as a function of log(a) = it/ln(10) for models with mo = 0.5 
eV and different values of mi and A. Bottom panel: Neutrino 
mass variation for the same parameters as in the top panel. 



transition is still on-going today, and the present epoch 
roughly coincides with the maximum variation. 

Although the functional form of T, eq. (fl~9]l . seems com- 
plicated, one should note that it is one of the simplest 
forms satisfying our requirements with a minimal num- 
ber of parameters. An example that could look simpler, 
but that for practical purposes is not, would be to assume 
that the two plateaus are linked together by a straight 
line. In this case, we would need a parameterization of 
the form 



mi 

m„ = ^ jtiq -f (mi — mo) 
m 



«NR — Wend 



U < 1t NR , 

!*NR < U < U cn d , 
U > Mend 



where u on d corresponds to the chosen redshift in which 
the transition stops. Notice that in this case not only 
we still have three parameters to describe the mass vari- 
ation, but also the function is not smoooth. Moreover, 
the derivative of the mass with respect to u gives a top- 
hat-likc function which is discontinuous at both unr and 
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itond- In this sense, it seemed to us that equation (|T8[) 
would give us the best "price-to-earnings ratio" among 
the possibilities to use phenomcnologically motivated pa- 
rameterizations for the mass- varying neutrinos, although 
certainly there could be similar proposals equally viable, 
such as for instance the possibility of adapting for the 
mass variation the parameterization used for the dark 
energy equation of state in [42|, EH- There, the tran- 
sition between two constant values of the equation of 
state exhibits a tanh [T t (u — u t )] dependence, where T t 
is responsible for the duration of the transition and u t is 
related to its half-way point. 

In the rest of our analysis, we will use a couple of extra 
assumptions that need to be taken into account when go- 
ing through our results. First, we will consider that only 
one of the three neutrino species is interacting with the 
dark energy field, that is, only one of the mass eigenstates 
has a variable mass. The reason for this approximation 
is twofold: it is a simpler case (compared to the case 
with 3 varying-mass neutrinos), since instead of 6 extra 
parameters with respect to the case of constant mass, we 
have only 2, namely the early mass of the neutrino whose 
mass is varying, mi, and the velocity of the transition, 
related to A. 

Besides simplicity, the current choice is the only one 
allowed presently in the case in which neutrinos were 
heavier in the past. Indeed, we expect our stronger 
constraints to come from those scenarios, especially if 
the neutrino species behaves as a nonrelativistic com- 
ponent at the time of radiation-matter equality, given 
by 1 + z cq - 4.05 x 10 4 (n c0 h 2 + n b0 h 2 )/(l + 0.23iV eff ) 
(here the indexes c and b stand for cold dark matter 
and baryons, respectively, and N e g is the effective num- 
ber of relativistic neutrinos). Taking the three neutrino 
species to be nonrelativistic at equality would change 
significantly the value of z oq , contradicting CMB data 
(according to WMAP5, 1 + z cq = 3141+^ (68% C.L.) 
j44|). Instead, a single neutrino species is still marginally 
allowed to be non-relativistic at that time. 

To simplify the analysis, we also assumed that the dark 
energy field, when not interacting with the neutrinos, 
reached already the so-called scaling solution (see, e.g., 
[H and references therein), i.e., the dark energy equation 
of state Wrf, in cq. (|15|) is constant in the absence of inter- 
action. Notice however that when the neutrinos become 
non-relativistic the dark energy fluid receives the anal- 
ogous of the chameleon kicks we mentioned before, and 
the dark energy effective equation of state, eq. (fl"6]) , does 
vary for this period in a consistent way. 

The upper panel of Figure [2] shows how the density 
parameters of the different components of the universe 
evolve in time, in a typical MaVaNs model. The lower 
panel displays a comparison between mass-varying and 
constant mass models, in particular during the transition 
from mi to mo- As one would expect, far from the time 
of the transition, the densities evolve as they would do 
in the constant mass case. 




0.001 
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FIG. 2: (Color online) Top panel: Density parameters for the 
different components of the universe versus log(a) = uj ln(10) 
in a model with mi = 0.05 eV, mo = 0.2 eV, A = 10, and 
all the other parameters consistent with present data. The 
radiation curve include photons and two massless neutrino 
species, and matter stands for cold dark matter and baryons. 
The bump in the neutrino density close to log(a) = —0.5 is 
due to the increasing neutrino mass. Bottom panel: Density 
parameters for two different mass-varying neutrino models. 
The solid black curves show the density parameter variation 
for two distinct constant mass models, with masses m„ = 0.05 
eV and m v = 0.2 eV. The dashed (red) curve shows a model 
in which the mass varies from mi = 0.2 eV to mo = 0.05 eV, 
with A = 0.1, and the dotted (blue) line corresponds a model 
with mi = 0.05 eV to m = 0.2 eV, with A = 10. 



C. Perturbation equations 

The next step is to calculate the cosmological perturba- 
tion equations and their evolution using this parameteri- 
zation. We chose to work in the synchronous gauge, and 
our conventions follow the ones by Ma and Bertschingcr 
[45l | . In this case, the perturbed metric is given by 



ds 2 



-a 2 dr 2 + a 2 (6ij + hij) dx l dx-' 



(20) 



6 



In this gauge, the equation for the three-momentum of 
the neutrinos reads 12 511 



dq 



1 • 2 m 2 dp A dx l 
._ qhm ^ a _f3__ 



where, as in equation ([¥]), we define 

dhim„ dh\m„ ( dp a 



(3(a) 



d\na 



(21) 



(22) 



Since the neutrino phase space distribution [45|] can be 
written as / (x\q 7 nj,r) = f"(q) [l + f (x\ q, n h r)] , 
one can show that the first order Boltzmann equation for 
a massive neutrino species, after Fourier transformation, 
is given by 0, HH] 



^ -I,- I \ T 

— + i- n-k * 

or e 



fj — (k • fi 



2 h + Qfi \ din/ 



dlnq 



. a qk a 2 mldlnf° 

(n- k) — — bp^ , 

e q d In q 



(23) 



where r\ and /i are the synchronous potentials in the 
Fourier space. Notice that the perturbed neutrino en- 
ergy density and pressure are also going to be modified 
due to the interaction, and are written as 



This extra term comes from the fact that the comoving 
energy e depends on the dark energy density, leading to 
an extra-term which is proportional to (3. 

Moreover, if we expand the perturbation \& (k, q, n, r) 
in a Legendre series [45j , the neutrino hierarchy equations 
will read, 



\j/ 2 



where 



qk 



qk 
3e 



hd In/ 
6 dlnq 

2* 2 ) + k , 



qk 
{2£+l)e 



- 3* 3 ) - 



5 11 



1 qka 2 m 2 dlnf° 
6 e q z dmq 



(26) 



din/ 
dlnq 



+ U 



(27) 



(>Pu 

3Sp u 



d 3 q 



(2tt) 3 
d 3 q 



/°(e* + ^ 2 



(27T) 



-f° 

3 •' 



(24) 



0^** ) -(25) 



For the dark e nerg y, we use the "fluid approach" [46j 
(see also [13, HH, |49|). so that the density and velocity 
perturbations are given by, 



3H{w A ~ cj) (84 + 



(l + ^) ^ + | - ^ /3p (l - 3t£)& + (3p A {l - 3w v ) 



(28) 



H{1 - 3cl)+/3p u (l - 3w v )H(l - 3w A ) 
l + Pp v (l-Zw v ) 



c^-/3(l - 3w„) 



Wcf, 



-p<j>Sc/> — P^d,), 



(29) 



where the dark energy anisotropic stress is assumed to 
be zero [5(J, and the sound speed is defined in the 
frame comoving with the dark energy fluid [5l|. So, in 
the synchronous gauge, the quantity c? = Sp^/Sp^ is 
related to through 



*'5, 



P4> Of 
Pa k 2 



f>4> 64, 



In addition, from eqs. ([15]) and (|22|) , we have that 

P0 _ -3H(l + w A ) 
p A 1 + 0p„(l — 3w u ) 



(30) 



(31) 



IV. RESULTS AND DISCUSSION 



A. Numerical approach 



Equipped with the background and perturbation equa- 
tions, we can study this scenario by modifying the numer- 
ical packages that evaluate the CMB anisotropics and the 
matter power spectrum. In particular, we modified the 
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TABLE I: Assumed ranges for the MaVaNs parameters 

Parameter Range 

— 1 < W4, < —0.5 

mo < mo/eV < 5 

A -4 < log A < 2 

fj, -6 < logO + ) < 

-6 < log(^t_) < 



CAMB code 1 [HlKbased on CMBFast 2 routines. We 
use CosmoMC 3 |54| in order to sample the parameter 
space of our model with a Markov Chain Monte Carlo 
(MCMC) technique. 

We assume a flat universe, with a constant equation 
of state dark energy fluid, cold dark matter, 2 species 
of massless neutrinos plus a massive one, and ten free 
parameters. Six of them are the standard ACDM pa- 
rameters, namely, the physical baryon density ftboh 2 , the 
physical cold dark matter density fl c oh 2 , the dimension- 
less Hubble constant h, the optical depth to reionization 
^roion, the amplitude (A s ) and spectral index (n s ) of pri- 
mordial density fluctuations. In addition, we vary the 
constant dark energy equation of state parameter w<f> and 
the three parameters accounting for the neutrino mass: 
the present mass tuq , the logarithm of the parameter A 
related to the duration of the transition, and the loga- 
rithm of the ratio of the modulus of the mass difference 
over the current mass, log p, where we define 

= \mx-mo\ j M+ = ^ - 1 ,m 1 >m , 
m \n_ = l-^ ) m 1 <m . 

All these parameters take implicit fiat priors in the re- 
gions in which they are allowed to vary (see Tabic [T| . 

Concerning the last parameter, notice that we choose 
to divide the parameter space between two regions: one 
in which the mass is decreasing over time (p+) and one 
in which it is increasing (p-). We chose to make this 
separation because the impact on cosmological observ- 
ables is different in each regime, as we will discuss later, 
and by analyzing this regions separately we can gain a 
better insight of the physics driving the constraints in 
each one of them. Moreover, we do not allow for models 
with < — 1, since we are only considering scalar field 
models with standard kinetic terms. 

For given values of all these parameters, our modified 
version of CAMB first integrates the background equa- 
tions backward in time, in order to find the initial value 
of leading to the correct dark energy density today. 
This problem does not always admit a solution leading 



1 |http : Z/canb. inf o/| 

2 |http : //cf a-www . harvard . edu/~mz aldarr /CMBFAST/cmbf ast . html 

3 http : //cosmologist . inf o/cosmomc/ 
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FIG. 3: (Color online:) Marginalised ID probability distri- 
bution in the increasing mass case mi < mo, for the neutrino 
/ dark energy parameters: mo, log 10 [/x_] (top panels), w,/,, 
and log A (bottom panels). 



to well-behaved perturbations: the dark energy perturba- 
tion equations lp?8|) , (|2l))) become singular whenever one 
of the two quantities, p^ or [1 + /3p u (l — 3w v )], appearing 
in the denominators vanishes. As wc shall see later, in 
the case in which the neutrino mass decreases, the back- 
ground evolution is compatible with cases in which the 
dark energy density crosses zero, while the second term 
can never vanish. We exclude singular models by stop- 
ping the execution of CAMB whenever p^ < 0, and giv- 
ing a negligible probability to these models in CosmoMC. 
The physical interpretation of these pathological models 
will be explained in the next sections. For other models, 
CAMB integrates the full perturbation equations, and 
passes the CMB and matter power spectra to CosmoMC 
for comparison with the data. 

We constrain this scenario using CMB data (from 
WMAP 5yr @ i|, VSA CBI [H and 

ACBAR [58|): matter power spectrum from large scale 
structure (LSS) data (2dFGRS jM and SDSS 001); su- 
pernovae la (SN) data from [6lJ, and the HST Key 
project measurements of the Hubble constant (62l| 4 . 

Once the posterior probability of all ten parameters 
has been obtained, we can marginalize over all but one 
or two of them, to obtain one- or two-dimensional proba- 
bility distributions. We verified that the confidence lim- 
its on the usual six parameters do not differ significantly 
from what is obtained in the "vanilla model" [44j, and 



4 While this work was being finished, the SHOES (Supernova, HO, 
for the Equation of State) Team 1631 reduced the uncertainty on 
the Hubble constant by more than a factor 2 with respect to the 
value obtained by the HST Key Project, finding H = 74.2 ± 3.6 
km s -1 Mpc - 1 . However, since we are taking a flat prior on Hq, 
and our best fit value for Ho is contained in their la region, we 
do not expect our results to be strongly affected by their results. 
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TABLE II: Results for increasing and decreasing neutrino mass, using WMAP 5yr + small scale CMB + LSS + SN + HST 
data. 

(+)Region 95% (68%) C.L. (-)Rcgion 95% (68%) C.L. 

w,/, < -0.85 (< -0.91) < -0.87 (< -0.93) 

m (eV) < 0.28 (< 0.10) < 0.43 (< 0.21) 

log At+ < -2.7 (< -4.5) 

log At. < -1.3 (< -3.1) 

log A [-3.84; 0.53] ([-2.20; 0.05]) [-0.13; 4] ([0.56; 4]) 



therefore we only provide the results for the extra neu- 
trino and dark energy parameters (Figures [JJ [H [U [31 and 
Table HU). 



B. Increasing neutrino mass 

In this model, the background evolution of the dark 
energy component obeys to equation (|15[) . which reads 
after division by p^: 

& = -3H(l + W(f> )-^^^n(l-3w u ) (32) 
p<j> du 

= -Td-r, 

where the two positive quantities Td and I\ represent re- 
spectively the dilution rate and interaction rate of the 
dark energy density. For any parameter choice, p^ can 
only decrease with time, so that the integration of the 
dark energy background equation backward in time al- 
ways find well-behaved solutions with positive values of 
Pj,. Moreover, the quantity [1 + (3p v {l — 3w v )\ appear- 
ing in the denominator of the dark energy perturbation 
equations is equal to the contribution of the dilution rate 
to the total energy loss rate, Td/ (Td + r,). This quantity 
is by construction greater than zero, and the dark energy 
equations cannot become singular. However, when the 
the interaction rate becomes very large with respect to 
the dilution rate, this denominator can become arbitrar- 
ily close to zero. Then, the dark energy perturbations 
can be enhanced considerably, distorting the observable 
spectra and conflicting the data. Actually, this amplifi- 
cation mechanism is well-known and was studied by vari- 
ous authors [2(1 HH, [65[ . It was found to affect the largest 
wavelengths first, and is usually refered as the large scale 
instability of coupled dark energy models. The condition 
for avoiding this instability can be thought to be roughly 
of the form 



Ti < AT, 



where A is some number depending on the cosmologi- 
cal parameters and on the data set (since a given data 
set tells how constrained is the large scale instability, i.e. 
how small can be the denominator [1 + @p v (l — Zw v )\, 
i.e. how small should the interaction rate remain with 



respect to the dilution rate). The perturbations are am- 
plified when the denominator is much smaller than one, 
so A should be a number much greater than one. In- 
tuitively, the condition (j39|) will lead to the rejection of 
models with small values of (w,/,, A) and large values of 
Indeed, the interaction rate is too large when the 
mass variation is significant (large p-) and rapid (small 
A). The dilution rate is too small when w$ is small (close 
to the cosmological constant limit). Because of that, it 
seems that when the dark energy equation of state is al- 
lowed to vary one can obtain a larger number of viable 
models if w<j, > —0.8 early on in the cosmological evolu- 
tion [is3.[67j. 

We ran CosmoMC with our full data set in order to see 
how much this mass-varying scenario can depart from a 
standard cosmological model with a fixed dark energy 
equation of state and massive neutrinos. In our param- 
eter basis, this standard model corresponds to the limit 
log/Lt_ — ► — oo, with whatever value of log A. The ob- 
servational signature of a neutrino mass variation during 
dark energy or matter domination is encoded in well- 
known effects, such as: (i) a modification of the small- 
scale matter power spectrum [due to a different free- 
streaming history], or (ii) a change in the time of mat- 
ter/radiation equality [due to a different correspondence 
between the values of (u>b, ui m , io v ) today and the actual 
matter density at the time of equality]. On top of that, 
the neutrino and dark energy perturbations can approach 
the regime of large-scale instability discussed above. 

Our final results - namely, the marginalized ID and 
2D parameter probabilities - are shown in figures [3J and 
[U The shape of the contours in (log/z_,logA) space is 
easily understandable with analytic approximations. The 
necessary condition (f33|) for avoiding the large-scale in- 
stability reads in terms of our model parameters 



i + A(i + r) 



< A 



1 



(i-r)(i-/) 



3O (l 



Wcf,) 



(33) 

where we expressed the mass variation as 



fi„(l - 2>w v ) 
(34) 



dlnm^ 



P- 



1 + A 



(i-r)(i-/) . (35) 



du \1 — P—T J V ^ 
Two limits can be clearly seen from this equations. For 
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FIG. 4: (Color online) Marginalised 2D probability distribution in the increasing mass case mi < mo. 



A <g; 1 (fast transitions), the upper limit on reads 



< AA 



(i-r)(i-/) 



3f^(l + w^) 
f2„(l - 3w v ) 



(36) 



This corresponds to the diagonal limit in the lower half of 
the right upper panel of figure [U In fact, the appearance 
of the large-scale instability is seen in models localized at 
the edge of the allowed region, as shown in figure El 

In the opposite case of a very slow transition, A» 1, 
it is clear from cq. (|3"4"|) that the limit on /x_ should be 
independent on A, 



t_ <A 



(i-r)(i-/) 



3^0(1 + W<j,) 

£!„(! - 3w„) 



(37) 



This limit corresponds to the almost vertical cut in the 
upper part of the plane (log p,- , log A) (upper right panel, 

fig- ED- 

These conditions are easier to satisfy when at the time 
of the transition, $7^(1 + w<p) is large. So, in order to 
avoid the instability, large values of are preferred. 
However, it is well-known that cosmological observables 
(luminosity distance relation, CMB and LSS power spec- 
tra) better fit the data for w close to —1 (cosmological 
constant limit). In the present model, the role of the 
large-scale instability is to push the best-fit value from 
-1 to -0.96, but Wcf, — — 1 is still allowed at the 68% C.L. 

The main result of this section is that the variation of 
the neutrino mass is bounded to be small, not so much 
because of the constraining power of large-scale structure 
observations in the regime where neutrino free-streaming 
is important (i.e., small scales), but by CMB and LSS 
data on the largest scales, which provide limits on the 
possible instability in DE and neutrino perturbations. 

Indeed, for the allowed models, the mass variation 
could be at most of order 10% for masses around 0.05 
cV, and less than 1% for masses larger than 0.3 eV: this 
is undetectable with small scale clustering data, showing 
that the limit really comes from large scales. 



With those results, we conclude that there is no ev- 
idence for a neutrino mass variation coming from the 
present data. In fact, as for most cosmological data anal- 
yses, the concordance ACDM model remains one of the 
best fits to the data, lying within the 68% interval of this 
analysis. 

Nonetheless, better constraints will possibly be ob- 
tained with forthcoming data, especially the ones that 
probe patches of the cosmological "desert" between z ~ 
1100 and z ~ 1, like CMB weak lensing [68j|, and/or 
cross-correlations of different pieces of data, like CMB 
and galaxy-density maps (6f|. We can estimate, for in- 
stance, what is the favored redshift range for the neu- 
trino mass variation according to our results. Taking 
mo = 0.1 eV and the mean likelihood values for log A 
and log [mi /mo], one can see that the bulk of the mass 
variation takes place around z ~ 20, a redshift that pos- 
sibly will be probed by future tomographic probes like 
weak len sing I70l . 17111 and especially 21 cm absorption 
lines [z3, [73|, 1 74 |75| . Those will help not only to dis- 
entangle some degeneracies in the parameter space, but 
will also allow for direct probes of the neutrino mass in 
different redshift slices. 



C. Decreasing neutrino mass 

In this case, the evolution rate of the dark energy den- 
sity is still given by equation (|32[) but with an opposite 
sign for the interaction rate: in can be summarized as 



— = —1 n 



(38) 



with and T, both positive. In principle, the interac- 
tion rate could overcome the dilution rate, leading to an 
increase of p^. Hence, the integration of the dark energy 
evolution equation backward in time can lead to nega- 
tive values of ps, and the prior pj, > implemented in 
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1 k (h Mpc" 1 ) 

FIG. 5: (Color on-line) CMB anisotropies and matter power spectra for some mass varying models with increasing mass, 
showing the development of the large scale instability. The cosmological parameters are set to our best fit values, except for 
the ones shown in the plot. The data points in the CMB spectrum correspond to the binned WMAP 5yr data. 



our CAMB version is relevant. Still, the denominator 
[1 + (3p v {l — 3uv)] can never vanish since it is equal to 

r d /(r d -ri). 

Well before before the transition, the interaction rate 
is negligible and p^ is always negative. We conclude that 
(3 = dh\m v / dprf, starts from small positive values and 
increases. If the condition 

r< < r d (39) 
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FIG. 6: (Color online:) Marginalised ID probability distribu- 
tion (red/solid lines) for the decreasing mass case mi > mo, 
for neutrino / dark energy parameters: mo, log[/x+] (top pan- 
els), w$, and log A (bottom panels). 



is violated during the transition, p$ will cross zero and 
become positive. This corresponds to (3 growing from 
zero to +00, and from —00 to some finite negative value. 
After Ti/Td has reached its maximum, /3 undergoes the 
opposite evolution. Reaching p§ = is only possible if 
Prf, has a non-monotonic evolution, i.e. if (|39p is violated. 
However, the perturbations diverge even before reaching 
this singular point: when (3 tends to infinity, it is clear 
from eq. (|26[) that the neutrino perturbation derivatives 
become arbitrarily large. We conclude that in this model, 
the condition (|39|) is a necessary condition for avoiding 
instabilities, but not a sufficient condition: the data is 
expected to put a limit on the largest possible value 1 of 
/3, which will always be reached before p^ changes sign, 
i.e. before the inequality (|39[) is saturated. Hence, the 
condition for avoiding the instability is intuitively of the 
form of (|33p . but now with A being a number smaller 
than one. 

We then ran CosmoMC with the full data set and ob- 
tained the marginalized ID and 2D parameter probabil- 
ities shown in figures [5] and [7J The major differences 
with respect to the increasing mass case are: a stronger 
bound on too, a much stronger bound on ^(_, and the 
fact that large values of A are now excluded. This can 
be understood as follows. In order to avoid instabilites, 
it is necessary to satisfy the inequalities (|57|) . but 

with a much smaller value of A than in the increasing 
mass case; hence, the contours should look qualitatively 
similar to those obtained previously, but with stronger 
bounds. This turns out to be the case, although in ad- 
dition, large A values are now excluded. Looking at the 
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FIG. 7: (Color online) Marginalised 2D probability distribution for decreasing mass, mi > mo. 



mass variation for large A in figure [TJ we see that in 
this limit the energy transfer takes place essentially at 
low rcdhsift. Hence, the interaction rate is large close to 
z = 0. In many models, this leads to positive values of 
Ptf, at the present time, to a non-monotonic behavior of 
the dark energy density, and to diverging perturbations. 
This can only be avoided when w is large with respect 
to -1, i.e. when the dilution rate is enhanced. Hence, 
in this model, the need to avoid diverging perturbations 
imposes a strong parameter correlation between w and 
A. However, values of w greater than -0.8 are not com- 
patible with the supernovae, CMB and LSS data set; this 
slices out all models with large A. 

The fact that the bound on mo is stronger in the de- 
creasing mass case is also easily understandable: for the 
same value of the mass difference /i± = \m\ — mo|/mo, a 
given m,Q corresponds to a larger mass mi in the de- 
creasing mass case. It is well-known that CMB and 
LSS data constrain the neutrinos mass through its back- 
ground effect, i.e. through its impact on the time of mat- 
ter/radiation equality for a given dark matter abundance 
today. The impact is greater when mi is larger, i.e. in 
the decreasing mass case; therefore, the bounds on mo 
are stronger. 



V. CONCLUDING REMARKS 

In this work we analysed some mass-varying neutrino 
scenarios in a nearly model independent way, using a gen- 
eral and well-behaved parameterization for the neutrino 
mass, including variations in the dark energy density in 
a self-consistent way, and taking neutrino/dark energy 
perturbations into account. 

Our results for the background, CMB anisotropics, and 
matter power spectra arc in agreement with previous 
analyses of particular scalar field models, showing that 
the results obtained with this parameterization are ro- 



bust and encompass the main features of the MaVaNs 
scenario. 

Moreover, a comparison with cosmological data shows 
that only small mass variations are allowed, and that Ma- 
VaNs scenario are mildly disfavored with respect to the 
constant mass case, especially when neutrinos become 
lighter as the universe expands. In both cases, neutrinos 
can change significantly the evolution of the dark energy 
density, leading to instabilities in the dark energy and/or 
neutrino perturbations when the transfer of energy be- 
tween the two components per unit of time is too large. 
These instabilities can only be avoided when the mass 
varies by a very small amount, especially in the case of 
a decreasing neutrino mass. Even in the case of increas- 
ing mass, constraining better the model with forthcoming 
data will be a difficult task, since it mimics a massless 
neutrino scenario for most of the cosmological time. 

One should keep in mind that our analysis assumes a 
constant equation of state for dark energy and a mono- 
tonic behavior for the mass variation. Even though those 
features are present in most of the simplest possible mod- 
els, more complicated models surely can evade the con- 
straints we obtained in our analysis. 

Finally, those constraints will improve with forthcom- 
ing tomographic data. If any of the future probes indicate 
a mismatch in the values of the neutrino mass at differ- 
ent redshifts, we could arguably have a case made for the 
mass- varying models. 
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